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INTERNAL PRESSURES OF A SPACECRAFT OR OTHER 
SYSTEM OF COMPARTMENTS, CONNECTED IN VARIOUS 
WAYS AND INCLUDING OUTGASSING MATERIALS, IN A 
TIME-VARYING PRESSURE ENVIRONMENT 

John J. Scialdone 
Test and Evaluation Division 

SUMMARY 

A method of calculating the pressure profile for each com- 
partment of a multi-compartment system such as a spacecraft is 
presented. The method is applicable to volumes connected in series 
or in parallel or both by passages of various sizes and geometries, 
where each volume may include several materials outgassing ac- 
cording to known or assumed characteristics. The entire system 
is exposed to a time-varying pressure environment. 

An appropriate set of differential equations with variable co- 
efficients is derived, which may be solved by computer to provide 
pressure profiles for each compartment. The profiles obtained 
are the result of viscous, intermediate, and molecular flow regimes 
experienced under appropriate conditions in the system of volumes 
and passages. Several computations which have been made are 
presented and compared to experimental results. Pressure pro- 
files for volumes with various outgassing materials and various 
orifice sizes have been calculated and are grouped according to 
the molecular flow regime time constants and the outgassing char- 
acteristics of the material. From these parametric plots, the 
quasi-steady pressure in a volume containing outgassing materials 
can be estimated. 
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INTERNAL PRESSURES OF A SPACECRAFT OR OTHER 
SYSTEM OF COMPARTMENTS, CONNECTED IN VARIOUS 
WAYS AND INCLUDING OUTGASSING MATERIALS, IN A 
TIME-VARYING PRESSURE ENVIRONMENT 


INTRODUCTION 

Many spacecraft failures have been attributed by inference, to improper 
pressure conditions (Refs. 1 and 2). For example, electrical voltage breakdown 
has occurred many times during spacecraft testing or during flight. In some 
cases it has occurred only during flight with systems which under testing ap- 
peared to be fre of this type of problem (Refs. 2 and 3). This breakdown occurs 
when the ambient gas attains its critical pressure, which is a function of the 
nature of the gas species, the type of electrodes used and the distance between 
them, number of charged particles present, and the magnitude of the electric 
fields. In addition to this difficult problem, some of the other malfunctions and 
problems dependent on pressure conditions are: 

• Loss of pressure in sealed components. 

• Loss of lubrication due to evaporation or sublimation. 

• deterioration of adhesives such as loss of plasticizer due to 
evaporation or sublimation. 

• Loss of convective heat transfer with resultant development of hot spots 
or cold spots. 

• Loss of dielectric resistance followed by corona discharge or arcing. 

• Fouling of terminals following condensation of outgassed material. 

• joss of system calibrat. on resulting from unavailability of convective 
heat transfer at the existing pressure. 

• Humidity, aerodynamic performance, radiation and ionization effects, 
electronic bombardment, etc. 

• Pressure difference across surfaces sufficient tc cause mechanical 
failure. 
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Many of the problems listed above can be eliminated or minimized by the 
use of pressurized containers, coatings, selective insulation, encapsulation, 
shielding, etc. (Refs. 3 and 4). Also, in many cases possible difficulties can be 
avoided expediently by employing selective commands for the experiment; e.g., 
operation of the system at a limited voltage can be commanded for unfavorable 
pressure conditions, with hill voltage operation at favorable pressure. Of 
course, it is hoped that those problems related to pressure conditions, as well as any 
other malfunction, will be discovered and eliminated during the thermal-vacuum 
testing of the spacecraft. These tests attempt to simulate the temperature and 
pressure encountered during launch and orbit flight. Chamber pressures below 
10 “ 6 torr are produced to establish the appropriate mode of heat transfer and 
hence provide simulation of the thermal environment. In testing, it is assumed 
that after long term vacuum exposure the pressure at any point in the space- 
craft will be comparable to that produced and measured in the chamber itself. 

But this cannot generally be verified, due to the absence of pressure instrumen- 
tation at inaccessible internal locations and due to the fact that many materials 
in the pressure range 10 ~ 5 to 10“ 9 torr expel adsorbed and absorbed gases. 

The latter as well as other gases produced by evaporation, distillation and de- 
composition, affect the pressure considerably. These ’’outgassing" conditions 
are likely to exist for a long time and set up, at certain locations of the space- 
craft, pressures higher than the chamber pressure (unbeknown to the experimenter). 
In addition, the transient pressures at the launch phase of the flight may produce 
pressure differentials depending on flow conductance, which cannot be reproduced 
during testing. Such differentials could establish conditions conducive to failures 
cr temporary malfunctions. 

In view of the above, it is the intent of this report to provide an analytical 
tool which can be used to estimate the variation with time ofpressures in a space- 
craft or other system of interconnected volumes. Knowledge of these pressures, 
which will permit undesirable pressure conditions to be identified, should be 
useful not only in connection with the qualification and reliability testing of the 
spacecraft, but also in determining the appropriate time for certain spacecraft 
operations. It is expected that the accuracy of the mathematical results will 
depend on the extent of knowledge of geometric leakage paths, container volumes, 
and material behavior. Many uses of this analytical tool will be indicated later. 

In general, the analytical tool must provide a description of the internal 
pressure vs time for each compartment of a system of compartments arranged 
in a most general way, and connected by passages of different types. The follow- 
ing data must be given or assumed in order to carry out the calculation: 

• External pressure vs time under flight conditions, or the chamber pressure 
during testing. 

• Type of ambient gas contained in the system. 
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• Temperatures of ambient gas and of internal surfaces. 

• Volume of compartments. 

• Surface areas of compartments and of components within the compart- 
ments. 

• Characteristics of the surface areas, i.e. the initial and final outgassing 
rates, and the outgassing time constant for each outgassing material. 

• Areas, dimensions and types of passages and leakage paths between 
compartments. 


ANALYTICAL MODEL 

Figure 1 shows a model of volumes and connecting passages chosen so as to 
represent the general case of a system of interconnected compartments. The 



p 0 f p 5 f p 4 f P| < p 2 f p 3 C 50 


Figure 1-Arrangement of Volumes and Passages Used as a Model for the 

Derivation of Flow Equations 
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flow equations derived for the model are intended to provide a flow description 
for any system. The model is arranged so as to provide both series and parallel 
flow through main passages, leakages between chambers, and outgassing from 
different surfaces. The prrows in the figure show the reference flow directions, 
corresponding to the relative magnitudes assumed for the chamber pressures. 
The acoial flow direction will be dictated by the actual pressure in each chamber, 
which will be a function ci the system parameters. Outgassing from surfaces is 
shown by arrows original ▼ at the separation walls. 

Pressure-Time History Equations 

Equating the change of mass of a gas stored in a volume V to the difference 
of the masses of gas entering and leaving that volume in time dt gives 


V dp = (m in -m out )dt . 


Since the change in density, according to the equation of state, is 



substitution yields 



= RT ( m in 


m 


out 


) • 


The flow rate, Q, may be defined in terms of a conductance, C, and the pressure 
differential, AP, as 


0 ^ CAP = RTm , 


so that the above isothermal equation of mass balance becomes: 
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V 


dP 

dt 


= C AP - C ft 

in 0 


AP . 


One can add to the right side of this equation the sources G. , representing the 
outgassing from the surfaces within volume i. This outgassing, as will be dis- 
cussed later, is independent of pressure but depends on temperature and on the 
nature of the outgassing materials in the chamber. The basic equation for the 
pressure history in a volume then becomes 



AP - C 0 AP + G . 


( 1 ) 


Applying this equation to each of the volume of Figure 1 in turn, for an assumed 
monotonic sequence of pressures 


P 0 - P 5 S P 4 5 P 1 - P 2 S P 3 ’ 

one obtains: • 

^1 - Y * + ^21 ) “ ^ l ) + ^31 (^3 " ^1 ^ 

V 1 

- < C 1S + L,s> ( p , - p 5 ) - L 14 ( p i - P 4 ) + 0, ] • 

^2 “ y ^^32 + ^ 32 ^ ^3 - ^2 ) ^21 + ^ 21 ^ ^2 “ ^25 ^2 “ ^ 5 ^ + ^ ’ 

2 

P 3 =h - < C 32 + L 32> < P 3 ' P 2> - L 3. ( P 3 ' P 1 > 

V 3 

- L 34 ( P 3 - P 4 )- L 33 ( P 3 - P 5 )+ G 3 )' 
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li «=-5' [L » (P >- P * >+L ‘< <P i- P < ) 

4 

-(C 4S +L 45 )( p , -p 5 ) +g 4 ], 

P 5 =V-[(C 15 + L 1S ) (P, -p 5 ) + L 25 (P 2 -P 5 ) 

V 5 

+ L 3 S (P3-P5)+(L 45 + C 4sHP4- P 5> 
-( C 50+ L 5 oH P 5-Po)+ G 5J- 


For a given P (t), the above equations can be generalized as follows to describe 
the rate of pressure change in the i th chamber of an N chamber arrangement: 


P. 



£ (Cij 4L ij )(P J -P,) , 

j =0 


( 2 ) 


where 


i = 1 , 2, 3, 


N; G. = 


- L 


-Gni’ 


1 J 


= C . and L. . = L 


j i 


n 


j i 


(The sum in 


L ^ 

n 
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is over all the different outgassing materials n contained in volume i.) The 
symmetry conditions for the conductances are valid in the molecular flow re- 
gime but not always in the viscous flow regime (Refs. & and 6). 

The sign of (P. - P ) determines the direction of flow between chamber i 
and chamber j. For the fictitious chamber arrangement shown in Figure 1, 
matrices could be used to provide the passage and leakage conductances. Any 
other chamber arrangement could also be described by equations (2) and the 
appropriate conductance matrices. 

A solution of equations (2), given the initial chamber pressures, the con- 
ductances, the environment pressure as a function of time, and the outgassing 
data, will provide the pressure -time history in each chamber of such an arrange- 
ment. The functional forms of the conductances, outgassing functions, and ex- 
ternal pressure forcing functions will be described next. 

Conductances 


The conductance of a passage is defined as the rate of flow of gas through 
the passage per unit of pressure difference producing the flow: 




( 3 ) 


The conductance is a function of the geometry of the passages, the tempera- 
ture, and the molecular weight of the gas flowing. It is also a function of the 
pressure and viscosity of the gas depending on the flow regime. 

The characteristics thatgovemthe flow regime of the fluid are the Reynolds 
number, R e = UpD/p, which relates the inertial force to the viscous force, and 
the Knuds en number, K n = A./D, which relates the mean free path of the gas to a 
characteristic dimension of the container (Refs. 5-7). 

For R e > 2300, the inertial force prevails over the viscous force, except 
for a narrow region adjacent to the walls, and the flow is turbulent. During 
evacuation, turbulent flow can occur near atmospheric pressure for a very short 
period. It is experienced when high capacity displacement pumps evacuate the 
system through improperly sized, long, small diameter passages (Refs. 6-9). 
This type of flow will not be included here because of its short duration, and 
(mainly) because of the difficulties it entails. 
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For R e < 2300, the viscous force prevails and the flow is viscous or 
laminar. 


The conductances employed here are appropriate for the three types of flow 
which occur at low pressure: viscous, slip-intermediate, and free molecular 
flow. The parameter which differentiates these regimes is the Knuds en number. 
This number, which can be also expressed as the ratio of the Mach number to 
the Reynold's number (Ref. 10), compares the mean free path of the molecules 
to a characteristic dimension of the container. It indicates whether the mole- 
cules are colliding with each other or with the wall of the container. The follow- 
ing criteria can be used to distinguish the flow regimes (Ref. 11): 

• For K n < 10“ 2 , the flow is governed by collisions of the molecules with 
each other and is viscous. 

• For 10" 2 < K n < 10, the flow is an intermediate type characterized by 
collisions of the molecules with each other and also with the container 
wall. 

• For K >10, the molecules are colliding with the confining wall and 
the flow is of the free molecular type. 

From the formulae for the conductances given below, it will be seen that the 
molecular conductances are independent of pressure. This is because the flow 
is a random motion of the molecules. The intermediate flow is a combination of 
the viscous and the molecular types of flow. The flow changes from completely 
viscous at high pressure, where a zero relative velocity between fluid and wall 
exists, to a transition flow with a relative velocity at the wall. At lower pres- 
sures, the flow becomes a free molecular flow. Conductances for free molecular 
flow are modified by an experimental coefficient to account for the transition 
flow region. 

Formulae for conductance of several geometrical arrangements are found in 
the literature (Refs. 5-3, 11-13). The list which follows includes the most com- 
mon conductances, taken, for convenience of presentation, fiom Ref. 8. These 
conductances, which are normally given for air at ambient temperature,, have 
been modified to make them applicable to other gases and temperatures. For 
geometric passages not listed, one can make use of the analogy between vacuum 
circuits and electrical circuits to obtain effective total conductances by combin- 
ing simple conductances in series or parallel, as follows: 


1 V* 1 

• Series combination arrangement: — = / . 

e L ' C i j 


(4) 
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• Parallel combination arrangement 


:C -=L c “ 


(5) 


1, Thin small aperture (Orifice) 


a. Viscous flow i.e. K < 10“ 2 

■ n 


For 1 > u > .52, 


C. } = 23.97^lJ /2 (1 -u) -1 u* 712 (1 - u ' 288 ) 1/2 A. 


( 6 ) 


For .52 > u > .1, 


/ T\ 1/2 , 

c i; = 6 * 27 (s) (i-ur 1 * 


1 J 


(7) 


For .1 > u > 0, 


Cij -6.27(1) 1/2 A. 


( 8 ) 


b. Molecular flow i.e. K > 10 

— — n 


© 1/2 

A 


(9) 


c. Intermediate Flow i.e. 10“ 2 < K <10 


C. . = 
» j 


(C n ) mol + (C n ) vise 


(10) 


-si 


S 


t 



2. Long pipe 


a* Viscous flow 

The conductance is inversely proportional to the viscosity, which 
in turn varies directly with the square root of the temperature 
(Refs. 11 and 14). The conductance of the air at 20°C has been 
modified by multiplying its coefficient by the ratio of the viscosities 
and by a temperature function, so that 


C i} = 3120 7] (T)~ l/2 


(P , +*i> D* 



( 11 ) 


where rj , the ratio of the viscosities, is given in Table 1. 
b. Molecular flow 




D 3 

L 


( 12 ) 


c. Intermediate flow 



(13) 


where the intermediate flow coefficient J is shown in Table 1. 
3. Short tubes 


a. Viscous flow 

One can use the series combination of orifice and pipe, i.e., 



1 

1 1 

(^ij)orifice (^ij)pipe 


(14) 
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Table 1 

Coefficients for Conductances (from Reference 8) 


H O 
2 


Mol. Wgt (grin) 

Vise. Air ~ 

” = Vise. Gas • 15 C 

2.02 

2.1 

4.0 

.93 

17.0 

1.86 

18.0 

1.90 

28.0 

1.04 

28.7 

1.0 

32.0 

.90 

39.9 

.82 

44.0 

1.24 


Intermediate Flow Factor J 


Rectangular Duct Flow 
Coefficients Y, a 


DP (torr-cm) 

2 

x 10' 2 

1 4 

x 10" 2 

6 

x 10" 2 

8 

x 10' 2 

1 

x 10" 1 

2 

x 10" 1 

4 

x 10 _1 

6 

x 10 ~ 1 

8 

1.0 

2.0 

4.0 

10.0 

x 10" 1 


1.1 

1.4 

1.7 
2.0 
2.3 

3.8 

6.9 

9.9 

13.0 

16.0 

31.0 

62.0 
153.0 


1.0 

.99 

.98 

.95 

.90 

.82 

.71 

.58 

.42 

.23 


1.108 


1.126 

1.151 

.198 

.297 

1.444 


Coaxial Tube Factor K 


/ D , 

















b. Molecular flow 


c u =3.8l( s ) 


T N l/2 


D 3 


L + 4 / 3D 


(15) 


c. Intermediate flow 


r> 


i j 


(J) 


( 16 ) 


4 . Rectangular Duct 

a. Viscous flow 

The conductance of air at 20°C, modified in the same manner as for 
the long pipe, y ields the following: 

= 4455 i7(T)-‘ 2 (Y)^ (- L y- i ), (17 > 

where the values of 17 and Y , a flow coefficient, are listed in Table 1 . 

b. Molecular fl_w 



= 9.68 (a) 



a 2 b 2 

L ( a + b ) 


(18) 


where the value of coefficient a is listed in Table 1 . 

c. Intermediate flow 
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5. Long annular tube 


a. Viscous flow 


t-i/2 

C.. = 3050 77 l — 


°1 _D 2 


(p; -uj) 

ln e D r /D 2 


(Pj 

2 


b. Molecular flow 


C. . 

i j 


3.81 



(K) 


(D t -D 2 ) 2 (D t + D 2 ) 
L 


(19) 


( 20 ) 


where the value of the coefficient K is indicated as a function of the 
ratio of the diameters in Table 1. 

c. Intermediate flow 





mol 


(J) • 


Initial Conditions - Boundary Conditions 


In addition to the geometric description of the system, one must specify the 
type of gas and the pressure and temperature conditions which exist at the start 
of evacuation of the system. It is also necessary to establish the pressures at 
which the different types of flow regimes exist (namely the viscous, intermediate, 
and molecular flow regimes). At these pressures, which are established by the 
passage dimensions and the mean free paths of the molecules, the conductances 
appropriate for each flow must be used for the integration of the flow equation. 
As indicated previously, the Knudsen number provides a criteria for the change- 
over (Ref. 11): For K n < 10 ~ 2 the flow is viscous; For K n > 1.0 the flow is 
molecular. 

The Knudsen number, with the substitution of the expression for the mean 
free path of the gas in terms of the pressure, temperature, and size of the mole- 
cule, becomes 



1 


BT 


( 21 ) 


K_ 



D V2 vn a 2 


2,331 x 10~ 2 °T 
D Y 2 7TCT 2 p DPc * 


where a is the molecular diameter of the gas ( a = 3.74 x 10 " 8 cm for air), n is 
the number of molecules per unit volume, and B is the Boltzmann Constant. 

The flow regimes can now be established in terms of the pressure P and a 
dimension D as follows: 


• Molecular flow: K >1.0, 

n 

DP < 2*331 x 10~ 20 T 

O ' 2 


(22) 


which for air at 25°C becomes 


DP < 5 x 10 3 ( torr - cm). 


(23) 


• Viscous flow: 


K n <10“ 2 , 


DP > 2.331 x 


10“ 18 T 


(24) 


which for air at 25°C is 


DP > . 5 ( torr - cm). 


(25) 


♦ 


In the above, P is an average pressure and D a characteristic dimension of the 
container. This dimension is, for example, the internal diameter of a pipe. In 
general, it can be taken as being the smallest dimension in a passage or as an 
equivalent diameter for a passage of area A, i.e., 


D = 



(26) 


14 


r 



The pressure of the pipe flow P, is taken to be the average of the inlet and exit 
pressures. 

Outgassing *: unctions 

Also affecting the ambient gas in the container, are gases which evolve 
from the container surfaces and the equipment inside the container under 
vacuum conditions. These gases are those which were physically and chemically 
sorbed by the surfaces, those which were in solution in the material, and those 
arising from evaporation and decomposition products. The exact mechanism of 
outgassing is still a matter of controversy, but theoretical and experimental de- 
terminations of outgassing rates have been reported in the literature (Refs. 15 
and 16). 

It is found from experience that the outgassing of materials may be 
characterized, in general, by the following statements: 

• The cutgassing depends on the nature of the materials; for a given material 
it is a function of the material’s surface treatment and of the atmospheric 
composition to which the material was exposed prior to exposure f o 
vacuum. 

• The outgassing is proportional to the surface area and increases rapidly 
with the temperature. It decreases slowly with time and apparently never 
becomes zero at ambient temperature. 

• At pressures lower than 10" 3 torr, the outgassing rate is independent of the 
pressure; i.e., the number of molecules expelled by the surface per unit 
time does not depend on the pressure above the outgassing surface. This 
behavior has not been established with complete certainty, but is valid as 

a first approximation. 

Outgassing rates of various materials at ambient temperature have been 
plotted as a function of time. Figure 2 shows degassing rates for some ma- 
terials. These outgassing fluxes have been obtained by exposing the materials 
to vacuum and comparing the rates with those provided by known conductances, 
or by reading the rates of pressure increase. 

Tabulation of outgassing properties of many materials is indicated in 
Tables 2 and 3. These tables, reproduced from Ref. lo for metallic, elastomeric, 
plastics, and ceramic materials at25°C,list the outgassing rates and their slopes 
at 1, 4, and 10 hours of exposure to vacuum. A knowledge of the outgassing rates 
of the materials at each instant is required for the solution of the pressure 



DEGASSING RATE (TORR-1 Aec/cm z ) 



Figure 2-Degassing Rate Per Unit Surface Area Versus Time for Some Materials 
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Table 2 

Outgassing Rate Constants for Metallic Materials 
(From Reference 15) 

Am >■ tret of sample (cm 1 ); 

5« » pumping speed for air at 25 *C (1/eec); 

Kk w air equivalent outgaaaing rata after A hr of pumping (Torrd»scc~ t -cm~' > ); 

•a « absolute value of elope of log-log graph of outgaaaing rata vs. time after A hr of pumping. 


Material 

Am 

s. 

10**4 


10*iC 4 

«4 

10 **„ 

«IS 

1. Stainless ateel EN 58B (polished, vapor degreased) 


1 



_ 

0.049 

1.6 

0.014 


2. Stainless steel 

65 

3.5 

1.75 

U 

0.44 

0.84 

0.21 

0.75 

3. Stainless steel 

— 

0.7 

0.9 

0.75 

0.3 

0.75 

0.2 

0.75 

4. Mild steel 

65 

3.5 

5.4 

1 

1.4 

1 

0.5 

1. 

5. Mild steel (slightly rusty) 

— 

1 

6.0 

3.1 

0.28 

1.5 

0.13 

1 

6. Steel, rusty 

65 

3.5 

44 

1.4 

5.5 

1.4 

1.6 

1.3 

7. Iron 

12 

— 

4 

1 

1 

1 

— 

— 

S. Nickel plated mild steel (polished, vapor degreased) 

— 

1 

5 

2 

0.036 

1.4 

0.01 

— 

Nickel plated steel 

65 

3.5 

2.8 

2 



— 

— 

10. Chrome plated mild steel (polished, vapor degreased) 

— 

1 

0.1 

1 

0.023 

0.93 

0.009 

— 

11. Aluminum, anodized 

65 

3.5 

— 


3.6 

0.8 

1.1 

2.3 

12. Aluminum, bright rolled (cleaned in Stergena) 

65 

3.5 

— 

— 

0.22 

1 

0.075 

1 

13. Aluminum 

12 

— 

15 

1 

3.7 



— 

14. Duraluminum 

— 

0.7 

1.7 

0.75 

0.6 

0.75 

0.35 

0.75 

IS. Aluminum spray coated mild steel 

65 

0.4 

0.6 

0.75 

0.2 

0.75 

I 0.1 

0.75 

16. Aluminum spray coated mild ateel (rusted) 

65 

0.4 

1.8 

0.65 

0.8 

0.75 

0.4 

0.75 

17, Copper (24 hr at 95 per cent humidity) 

100 

0.32 

0.2 

2 

— 

— 


1 — 

18. Copper 

12 

— 

23 

1 


. — 

— 

— 

19. Brass (24 hr at 95 per cent humidity) 

43 

0.32 

0.15 

2.5 

— 

— 

— 

— 

20. Brass, cast 

650 

3.5 

10 

1 

2.5 

1 

1.1 

0.75 

21. Brass, wavc*guide section 

50 

0.3 

4 

2 

0.3 

1.4 

0.1 

1.2 

22. Nickel 

12 

— 

6 

1 

1.5 

1 


* 

2T. Nickel 

12 

— 

10 

1 

2.5 

1 



24. Molybdenum 

12 

— 

7 

1 

1.7 

1 

_ 



25. Tantalum 

12 

— 

9 

1 



r -TT 

- 

26. Zirconium 

12 

— 

13 

1 





27. Tungsten 

12 

— 

2 

1 



... 



. 

78. Silver 

12 


6 

1 

— 

— 

*— 

— 


history equation. Equations for the outgassing as a function of time, G(t), have 
therefore been fitted to the data in the tables. Techniques described below were 
used to calculate the parameters of these € quations for several materials. 

For outgassing rates at temperatures other than 25°C, an estimate of the 
increase of the outgassing rates with the temperature can be obtained from the 
empirical relation proposed by Arrhenius for the increase of the rate of chemi- 
cal reaction with the temperature: 


Q = Q 0 exp 



where Q is the rate at temperature T, Q 0 a constant (determined from Q at 25°C), 
AH a constant representing the energy of the activation of the material, R the gas 
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Table 3 

Outgassing Hate Constants for Non-Metallic Materials 

(From reference 15) 


Am » area of sample^cm 1 ); 

Sa *=* pumping speed for air at 25 °C (1/sec); 

Kt *• air equivalent Out an Ming rate after 1 hr of pumping (Torr^l-sec^-cm-*) ; 

K a outgassing rate after 4 hr of pumping; 

a x ** absolute value of slope of log-log graph of outgassing rate v». time after l hr of pumping; 
a 4 absolute value of slope of log-log graph after 4 hr; 

» time at which a begins to increase rapidly above o t for given sample (hr); 
w* » square of sample thickness! xo m (mm). 


Material 

A m 

s « 

10 7 iC, 

«1 

107C, 

*4 

t m 

V>1 


A* Elastomers 

t . Natural gum rubber (32 Durometer) 

75 

0.5 

12 

0.5 

6.0 

0.5 

12 

9 


Natural gum rubber (32 Durometer) 

1 

0.5 

13 

0.5 

6.5 

0.5 

10 

9 


2. Natural white rubber J.1260 total outgassing 

100 

0.35 

11.5 

0.40 

6.6 

0.42 

S 

2.6 


3. Natural white rubber J.1260 non-condensable in liquid 

100 

0.35 

2 

0.5 

0.2 

1.6 

— 

2.6 


air 

4. Natural crcpc rubber (vulcanized with S 20 min) 

22.2 

0.7 

73 

0.7 

31 

0.65 

— 

• 


5. Natural crepe rubber (vulcanized with Te 20 min) 

22.2 

0.7 

61 

0.5 

28 

0.60 

— 

• 


6. Neoprene F-905 

65 

0.3 

50 

0.49 

25 

0.55 

>30 

16 


7. Neoprene 60° Shore 

— 

1 

80 

0.45 

43 

0.49 

> 7 

— 


8. Neoprene 

65 

3.5 

300 

0.4 

180 

0.4 

15 

— • 


9. Neoprene (24 hr at 95 per cent humidity) 

7.6 

0.32 

1400 

0.75 

480 

0.96 

> 5 

— 


10. Neoprene (outgassed-f 24 hr dry N f ) 

7.6 

0.32 

120 

0.5 

— 

— 

— 

— 


11. Neoprene GNA (V-4-2) (J hr 50 per cent rel. humidity) 

30 

2 

... 

— 

30 

— 

— 

— 


12. Neoprene (bell jar gasket, SO Durometer) 

68 

0.4 

30 

0.5 

14.5 

O.S 

— 

• 


13. Neoprene 1157 (sulfur free) 

65 

1.5 

54 

0.45 

30 

0.45 

>20 

10 


14. Red vacuum hose (24 hr at 95 per cent humidity) 

141 

0.32 

88 

0.6 

30 

0.8 

3 

• 


15. Red vacuum hose (outgassed+24 hr dry Ni) 

141 

0.32 

11 

0.5 

— 

— 

— 

• 


16. Pcrbunan 

65 

3.5 

35 

0.3 

22 

0.5 

5 

— 


17. Pcrbunan PD 651 (vulcanized 30 min) 

22.2 

0.7 

56 

0.4 

31 

0.5 

4 

• 


18. Pcrbunan DR-39 (low acrylonitrile) 

22.2 

0.7 

94 

0.65 

41 

0.56 

20 

• 


19. Perbunan (Bayer) 

12 

0.1 

20 

O.S 

. 10 

0.5 

— 

— 


20. Perbunan (Bayer) (tempered) 

12 

0.39 

4 

0.59 

1.2 

1.1 

— 

— 


21. Butyl GRI (V-3) (i hr 50 per cent relative humidity) 

30 

2 

— 

— 

20 

— 

— 

— 


22. Butyl (DR-41) (40% C black; vulcanized 12 min) 

22.2 

0.7 

15 

0.68 

4. 

0.64 

— 

• 


23. Butyl (IUJ. 12000) 

100 

0.35 

20 

0.64 

6 

0.42 

— 

2.6 


24. Butyl (BU. 12000) (non-condensable liquid air) 

100 

0.35 

1 

0.42 

0.55 

0.51 

— 

2.6 


25. Convascal 

65 

3.5 

14 

0.2 

9 

0.39 

— 

— 


26. Convascal (extruued | in. square) 

39 

0.4 

10 

0.5 

4.9 

0.6 

>50 

91* 


27. Convascal 

100 

0.4 

5.0 

0.46 

2.6 

0.4 

>30 

10 


28. Polyijocyanate rubber D-43 

22.2 

0.7 

280 

0.45 

127 

0.57 

10 

• 


29. Nygon 

65 

3.5 

130 

0.5 

65 

0.6 

5 

— 


30. Hycar H-50 (ASTM spec. SB-510) 

6.S 

1.5 

120 

0.42 

70 

0.45 

>50 

10 


31. Hycar H-50 (ASTM spec. SB-510) 

65 

1.5 

140 

0.48 

74 

0.5 

>96 

10 


32. Silicone rubber (Wacker R 60) 

12 

0.39 

70 

1.07 

17 

1.1 

— > 

— 


33. Silicone rubber (Wacker R 80) 

12 

0.39 

180 

1.0 

44 

1.2 

— 

— 


34. Silicone tubber (24 hr 95 per cent humidity) 

100 

0.32 

230 

0.65 

46 

1.3 

— 

— 


35. Silicone rubber (outgassed+24 hr dry N*) 

100 

0.32 

13 

0.5 

— 

i — 

— 

— 


36. Silastic 

— 

1 

25 

0.5 S 

6 

1.8 

— 

— 


37. Silicone rubber 

22.2 

i 0.7 

94 

0.75 

31 

0.8 

1 20 

• 


38. Silastic X-6145-C 

65 

1 1.5 

25 

1.0 1 

5.6 

1.07 

>96 

10 


39. Silastic 8-164 (red, 62 Durometer) 

65 

1.3 

12 

o.9 : 

3.7 

0.9 

>70 

10 


40; Silastic **.0 (white, cured 24 hr at 480 °F, 74 Durometer) 
41. Silastic 50 (white, 55 Durometer) 

65 

1.3 

28 

1.0 

6.0 

1.0 

>50 

10 


05 

1.5 

3t> 

1.0 

6.4 

1.0 

>50 

10 


42. Silastic 67-163 (red, 61 Durometer) 

65 

1.3 

19 

0.93 

5.4 

1.0 

>50 

10 





18 


I 


Table 3 (continued) 


Material 


B. Plastics 

1. Teflon, Dupont 

2 . Teflon, high temperature (24 hr at 95 per cent 

humidity) 

3. Teflon, high temperature (outgassed 4* 24 hr dry Nt) 

4. Teflon, Ccroc (24 hr at 95 per cent humidity) 

5. Teflon, Dupont 

6. PTFE 

7. Kel-F (Oak Ridge National Laboratory) 

8. Kel-F 270 (Resistoflex Corp.) 

9. Araldite D 

10. Araldite D 

11. Araldite B 

12. Araldite F 

13. Araldite, Type 1 (Ciba, cured 15 hr 266 °F) 

14. Epoxy resin 200 (24 hr 95 per cent humidity) 

15. Epoxy resin 200 (outgassed +24 hr dry N t ) 

16. Eppon, Shell Oil Co. 

17. Plexiglas (Alsthom) 

18. Plexiglas M222 (Rohm and Haas) 

19. Plexiglas (24 hr at 95 per cent humidity) 

20. Plexiglas (outgassed + 24 hr dry N,) 

21. Mylar (24 hr at 95 per cent humidity) 

22. Mvlar V-200 (24 hr at 95 per cent humidity) 

23. Mylar V-200 (outgassed + 24 hr dry Ng) 

24. High temperature Thermalon (24 hr at 95 per cent 

humidity) 

25. Polyvinylchloride (24 hr at 95 per cent humidity) 

26. Polyvinylchloride (outgassed + 24 hr dry Ng) 

27. Polyethylene (BASF) 

2S. Polyethylene (Dynamit-AG) 

29. Polvamid (Payer) 

30. Ultramid (BASF) 

31. Poly viny lea rhazol (BASF) 

32. Poly sty rol (BASF) 

33. Poly sty rol (Dynamit-AG) 

34. Polyurethane (Bayer) 

35. Tcxtolite 11564 (24 hr at 95 per cent humidity) 

36. Tcxtolite 11564 (outgassed + 24 hr dry N,) 

37. Nylon 

38. Pcramfil (General Electric Co.) 

39. Polyester (Zenith Plastic Co., fiber glaaa laminate) 

40. Polyester (Plastone, fiber glaaa laminate) 

41. Polyester (Perault) 

42. Polyester, Norsodyne 

43. Celluloid 

44. Polythene 

45. Methylmethacrylate 

46. Epicote 

C. Ceramic* 

1. Porcelain, glased 

2. Steatite 


BB 

K9 

10’JC, 

n 

lO’/C. 

•# 

tm 

Wm 

12 

0.23 

n 

0.7 

1.6 

0.7 

H 

rnrnmm 

27 

0.32 

El 

0.93 

6 


§gl 

— 

27 

0.32 

0.1 

2 




— 

42 

0.32 

35 

0.46 

19 

— 


— 

65 

1.5 

5 

0.68 

1.2 

1.2 

3 

2.6 

— 

1 

3 

0.45 

1.5 

ESI 


— 

65 

1.5 

0.4 

0.57 

0.17 


I 

10 

47 

0.4 

— 

— 

0.16 

Ell 


0.6 

650 

3.S 

80 

0.8 

22 

0.78 

1 

— 

30 

0.7 

19 

0.3 

12.5 

0.5 

1 

— 

30 

0.7 

18 

0.4 

9.2 

wm 


— 

30 

0.7 

15 

0.5 

7.3 



— 

62 

•0.4 

5.5 

1.3 

0.9 

1.3 

10 

0.3 

50 

0.32 

110 

0.6 

— 

— 

— 

— 

50 

0.32 

0.2 

1.6 

— 

— 

— 

— 

41 

0.4 

40 

1.2 

8.5 

1.2 

50 

40 

30 

0.7 

31 

0.4 

18 

0.4 

_ 

— 

12 

0.10 

19 

0.43 

10 

0.57 

_ 

— 

62 

0.32 

115 

0.65 

40 

0.75 

— 

— 

62 

0.32 

10 

05 

— 


— 

— - 

289 

0.32 

0.1 

4.3 

— 

— 

_ 

— 

370 

0.32 

23 

0.7S 

4 

« 

— 

— 

370 

0.32 

1.3 

1.33 

— 

— 

— 

— . 

32 

0.32 

13 

2.3 

0.1 

— 

— 

— 

.70 

0.32 

8.5 

1.0 

0.2 

_ 



— 

70 

0.32 

0.08 

3.2 

— 

— 

— 

— 

12 

0.10 

2.3 

0.5 

1.15 

o.s 

— 

— 

12 

0.23 

2.6 

0.5 

1.3 

0.5 

— 

— 

12 

0.12 

46 

0.5 

23 

0.5 

— 

— 

12 

0.1 

17 

0.5 

8.5 

0.5 

- - 

— 

12 

0.23 

16 

0.5 

8 

0.5 

— 

— 

12 

0.14 

6 

0.5 

3 

0.5 

_ 

— 

12 

0.14 

15 

I 0.5 

I 7.5 

0.5 

— 

— 

12 

0.85 

5 

0.5 

1 2.5 

0.5 

_ 

— 

61 

0.32 

55 

0.7 

18 

1.1 i 

_ 

— 

61 

0.32 

0.1 

3.3 


WSM 

— 

i — 

— 

1 

120 

0.5 

60 

KXB 

— 

1 — 

57- 

0.4 

400 

-1.8 

1 40 

1.2 

50 

I 23 

52 


23 

0.84 

7 

0.84 

50 

16 

59 


25 

| 0.84 

8 

0.81 

50 - 

40 

30 


34 

0.72 

13.5 

0.67 

40 

40 

430 

3 jfl 

16 

0.36 

10 

0.36 

100 


12 


86 

0.5 

45 

0.5 

1 — 

— 

80 

rffSm 

200 

| 1.6 

i 20 

1.6 

— • 

— 

174 

0.7 

1 42 

I 0.9 

14 

0.57 

— 

— 

82 

0.7 

25 

m 

12.5 

0.5 

— 


30 

0.7 

6.5 


3.0 

I 

mmm 

mmm 

30 

0.7 

0.9 

i 

0.24 



' 
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constant, and T the absolute temperature. The rates of many organic reactions 
are doubled for each 10°K rise in temperature. A rule of thumb frequently sug- 
gested is that a 10% change in absolute temperature will cause a 10:1 change in 
outgassing rate (Ref. 7). 

From the plots of Figure 2, the outgassing rates are seen to be represen- 
table by the following functions: 


* = % t_a 


(27) 


and 


q = q 0 e” t / T + q f , 


(28) 


In these expressions q Q is the initial outgassing rate, q f is the final out- 
gassing rate, a is the slope of the outgassing curve (27) and r is a time constant 
which represents, for the curve (28), the time required for the outgassing rate 
to drop to a value 



368 ^ + q f . 


It should be pointed out that an expression such as the above could be ob- 
tained experimentally for an entire package containing various materials. The 
equivalent outgassing so obtained could then be used in calculations involving 
the interaction of this package with others. The outgassing can be programmed 
as a mathematical function or provided in tabular or graphical form for the use 
to be made of it here. 

*• 

With the outgassing rates at 1 hour, 4 hours and 10 hours provided by the 
tables, the parameters of the exponential equation (28) representing the out- 
gassing can be obtained as follows. Using the system of equations 



Q 4 = Ae*‘ 4/C + B, 


Q 10 = Ae- 10 / C B| 


where the Q A are the known rates, and A, B, C are the parameters to be found, 
the following ratio is obtained: 


1 - ^4 ~ Qio _ e~ 4/c ( 1 - e~ 6/c ) _ 1 + e~ 3c = e - 3 /c + e -6/c (29) 

R " Q, -Q 4 e“ 1/c ( 1 - e” 3/c ) ’ e 3/c 

Letting x -1 = e ” 3/c and x -2 = e“ 6/c , this may be rewritten 


— = — , or x 2 -Rx-R = 0, 

R x 2 


which when solved for the positive root yields 


.. . r 3/c _R+ y R 2 + 4R 




from which the parameter 1/C is found: 


1 In x 
C~ 3 



R + ^R 2 + 4R 

2 


(30) 


Si 
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With C known, the parameter A can then be obtained as 


and finally 


A = 


Qi 



( e" 4/c - e“ 1/c ) * 


B - Q 10 - A e' i0/c > 


( 31 ) 


(32) 


The tables provide outgassing rates at only two different times for some 
materials. However, the slopes of the outgassing, a , are provided at these 
points. In these cases, to fit the data to the exponential equation, one can use 
the average of the slopes to establish 


l'C 


a. 



(33) 


The parameters A and B are then obtained as above, by using equations like (31) 
and (32). 

Table 4 presents the values of A, B, and C calculated for some of the ma- 
terials listed in Tables 2 and 3. 

If the tables provide only one rate of outgassing and a corresponding slope, 
the outgassing equation to be used is Q = At “ a . The coefficient A can be obtained 
directly by substitution of the given values, i.e., 


a* 


In general then 


Q = (Q, 




-a 


(34) 


► 
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Table 4 

Coefficients for the Outgassing Equations 


A 

a 

2 . 093507046-0 7 

8 * 2.04534053F-08 C * 1 . 594460276 00 

stainless 



check 

CHECK 

1. 75000000E-07 1 . 75000000E -0 7 

^4OUQD63OE-O0 . 4.400000006-08 



CHECK 

Z.lOliOOOOOE-OS 2. 1000000 06- 00 


A 

s 

i*3Jboist«e*o/ 

CHECK 

8 * 1 • 978372 7?F* 08 C * 1.550340896 00 

9.noaoncOQ£-00 O^DUOOiifliiaE-OB 

stainless 



CHECK 
CHECK. _ 

3.0000nOOOE-08 3, 0000000 HE *08 

2. OOOOOOOOE -06 2.00000000E-08 


A 

s 


k s 4«645b6G!iAE-Jl4 C * l.Bpib69l«E OP 

KiLUS 1 EEL 



CHECK 

CHECK 

5.40000000E-07 5.400000006-0/ 

1. 4 0 0 0 0 0 0 0E- 0 7 1.400000006-0/ 




CHECK 

8 . 000000008*08 5 , 000 0000 0E- 00 


A 

s 

4 .420837796 03 
CHECK 

8 * 1 • 303681608*00 C * 1.814617986*01 

6-ODOOOOOOE 00 fc.ooooonunE _JDJJ 

MILOS 1 66l 



check 
LtiEZIL 

2. 800000006* *! 8 2.0C38516OE-O0 



a 

9 t 3/5fl7lU>E-Q6 

k :: 1^56*1 y2J9£*>A7 , „ , . C s 00 

5iEctiJSlY 



CHECK 

CHECK 

4.40000000E*U6 4,4000000"E*06 

0.0 0 0 Q (1 0 0 Q£ - U 7. . 5.800000006-0/ 



CHECK 

1.600000006-07 1.600000006*0/ 


A 

» 

2-OUb29l5lE*06 

CHECK 

8 * 9.999998898*08 C « 5 .25967703E‘O1 

A.IUHLOIULOOE-07 4. 000000006*0/ 

* IRON 



CHECK 

CHECK 

1.O1000000E-07 1.010000006-0/ 

1.0D300U006-S7 Ju.0 Q_0 J).0 J) QJ)£ - 0 7 



* 


k s 9.999^9^17^-10 C = b . 

NiCKEwPQL . 



■'js« 

.CHECK 

8 . 0 000000 0E* 0 7 8 . OOOOOOOOE- 0 / 

3.600000006-09 3.600000016-OV 



CHECK 

1. OOOOOOOOE-09 1.00000001E-09 


7^ 

* 

1 • 6929505<6-08 
-CHECK 

8 * 8- 645897039*10 C * 1-62101720E 00 

j. ^ 0.0 ojo nooofc-oe i.oooooouoe-ob 

CHRUMtPLATESIEEL 



CHECK 
CHECK ... 

2.30000000P-09 2 . 300 0000 0E- 09 

9 JLaJ3.DOOO.ttE- 10 9. OOOOOOOOE-IO 


A 
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where Q n , t n and a are the values given in the tables. In particular, if the given 
value of Q is that corresponding to one hour, then A - Q . 

n n 

Ex ternal Pressure Ig (t) 

The pressure of the environment must be known as a function of time in order 
to solve the system of differential equations. Two vacuum environment profiles 
will be encountered by a spacecraft - the pressure versus altitude during launch 
and orbiting, and the pressure produced by pumping devices in a test chamber. 

The launch and orbiting pressure versus time profile for a spacecraft can 
be obtained by combining the pressures versus altitude as found in tables with 
the altitude versus time of the spacecraft (Ref. 17). For example, an equation 
for the external pressure can be written using simplifying assumptions, for a 
satellite having a period T # in an elliptical geocentric orbit. Let P max and 
P . ( be the environmental pressures at perigee and apogee respectively, P av 
the average pressure during orbiting and t Q the time from launch to reach orbit 
altitude. Then one may write simply: For 


while for 


where 


0 < t v t 


o ’ 



t > t 


o * 




max 


min 


max 


-p 


mm 


s i n oj t , 


(35) 


(36) 


^ = 27r/T # , and P atm = atmospheric pressure at launch . 


In the above equations, it has been assumed that the drop in pressure is 
linear with altitude during launch (which is approximately true for altitudes up 
to 125 km). For computing purposes, however, it will only be necessary to have 
pressure versus time in tabular form. 

The evacuation of a test chamber is a function of its volume and the pumping 
speed of the system. It is dependent on the temperature, humidity, and surface 
conditions of the chamber, its material, the gases present and many ether factors. 
It is thus quite possible that the results of one test will be somewhat different 
than another. The pressure versus time data for a particular chamber can be 
described analytically and used for the solution of the chamber-spacecraft 



system. As an example of this, Table 5 shows the pumpdown equations for an 
8' x 8' chamber, in part derived from the data shown in Figure 3. The data at 
the 55°C and * 10°C test chamber temperatures shown in this figure cover only 
pressures below 10“ 3 torr. The evacuation time from 760 torr down to 10 " 3 
torr is a small fraction of the total evacuation time and was not shown on this 
graph. The pressure for this interval was assumed to decay exponentially from 
time zero to .9 hr., followed by a straight line decay to 1 hour. This is in ac- 
cord with other experimental evidence. For the same time range, a straight line 
pressure drop can also be assumed. Calculations to fit these data and the above 
assumptions have provided the equations shown in Table 5 for the empty 8' x 8' 
test chamber at the two different chamber wall temperatures. 

Equations for chamber pressure profiles like these indicated here can be 
derived for other chambers. However, test data on the evacuation of any 
chamber can also be used directly in tabular form for the solution of the pres- 
sure equations. The latter method will generally be used here to calculate 
pressure profiles. 

Heat Transfer Considerations 


The pressure relation equation (2) has been derived under the assumption 
of isothermal flow conditions. This presupposes an average speed of evacuation 
and sufficient heat transfer from the container to the ^as to maintain a constant 
temperature. This condition is fact approached during evacuation of a test 
chamber or in actual flight, with the exception oi the very start of evacuation. 

At the start, there is a very rapid evacuation to represent the pressure drop 
experienced during rocket ascent. The flow is turbulent and viscous, and is 
provided by the rapid changing environment or by large capacity displacement 
pumps. The rapid evacuation produces a sudden drop in temperature, 4-5 °C for 
a large vacuum chamber, followed by a temperature recovery L v . ok to the ambient 
wall temperature. During this very short and rapid exhaus*, • .h i heat transfer 
from the wall is not sufficient to maintain the gas temperature constant and 
adiabatic expansion of the gas is approached. When the rate of expansion of the 
gas becomes moderate, and the rate of pressure drop is reduced, the heat input 
from the wall is sufficient to compensate for the heat loss caused by the expan- 
sion. An estimate of the temperature effect can be carried out as follows. 

The heat gain of the gas i i time dt is 


dy = h A i T -T) dt , 

Xg 8 V W g ' 7 


( 39 ) 
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Table 5 

Pressure Profile Equations for an 8* x 8’ Chamber 
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where A s is the wall surface, T w and T are respectively the wall and gas 
temperatures, and h is a heat transfer coefficient covering conductive, convec 
tive and radiative heat transfer and dependent on many variables such as flow 
regime, gas type, and initial temperature. 

The heat gain of the wall in the same period is 


dC, 


<V - T »)= C ( T .'-T„) 


(40) 


where C is the thermal capacitance of the container and T w ' is the container 
temperature before it changes to T w . Setting dQ =- dQ w , the heating; or cooling 
rate of the wall due to gain or loss of heat to the gas is found to be 


C ^= hA , (T.-V • (41) 


From this equation, the gas temperature can be obtained provided .the coefficient, 
h, the initial conditions, and the way in which the wall temperature varies with 
time are known. The solution would have to take into account heat gained or lost 
by the walls, including the gains due to solar radiation, aerodynamic heating, 
albedo and IR radiation, and power sources in the system together with the losses 
due to veradiation, conduction and convection from the walls. 

One of the difficulties resides with the heat transfer coefficient. This co- 
efficient changes radically from the viscous to the molecular flow regime and 
is related to the gas pressure by an "S" shaped curve. As shown in the kinetic 
theory of gases, the coefficient of heat conduction in the viscous flow regime 
is directly related to the gas viscosity and to the coefficient of heat capacity at 
constant volume (Refs. 10, 11, 14): 


k = i^UC v \ = i,C v . 


(42) 


Consequently, it varies directly with the temperature and is not affected by the 
pressure. Making appropriate substitutions in this equation and including the 
variation of specific heat with temperature leads to the Sutherland relation for 
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the coefficient of heat conductivity for various gases (Ref. 11). One could use 
this coefficient of heat transfer at pressures in which viscous flow occurs. In 
this range, heat transfer by radiation is a negligible part of the total heat trans- 
fer because the emissivity of the internal surfaces is low, the temperatures in- 
volved are low, and the coefficient of conductivity is sufficiently high. For pre- 
sure corresponding to the molecular flow regime, i.e., when K > 1, the heat 
transfer becomes orders of magnitude smaller than the viscous coefficient 
(Refs. 18 and 19). The heat transfer is then provided by radiation, and conduc- 
tion and convection are almost inexistent. Depending on the dimension of the 
container, the changeover from high to low conduction heat transfer occurs be- 
tween 1 torr and 10 ~ 5 or 10“ 6 torr. 

From the above, it appears that as first approximation one could assume a 
coefficient of heat transfer varying inversely as the gas density. 


--•O' 

and that the heat balance equation would be 


( V - T. ) = h 0 Q A ( T„ - T, ) dt. (44) 

In order to get the ratio of densities appearing in (44), one must estimate the 
gas temperature, which will change with the pressure. From the energy equa- 
tion, assuming no inflow and the absence of external heat or work, the change in 
the total internal energy of the gas is equal to minus the enthalpy carried out by 
the gas, i.e., 


E' - E = AH , 
or 

C v tm' T,' - m T, ] = - (m - »*) C„ T,. (45) 
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Rearranging, dividing by m' , T and C y and making the substitution C p / C v = 7 , 


T ' 

s m , 1 \ 

t (i-y) + ? • 

1 m 
8 


(46) 


Assuming the perfect gas law and that the process is isochoric, 


PV . , P'V 

m = - — and m , 

TR RT ' 

g s 


one has 


P T * 


m _ r "g 

m ' " p ' T 


g 


which substituted into (46) gives 


T ' 

g 


7 


f g 1 + [(7 - 1) ^7 

P 


(47) 


The ratio of masses in terms of pressure is 


m' 1 fp' 1N 1 

- = — — + (7 - 1) • 
m 7 |_p 


(48) 


and because p = m/V and p' = m' /V the following is obtained for the density 
ratio in terms of pressure ratio: 


=![£.'+ (7 -1)1. 
p y\j> J 


( 49 ) 
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This expression allows one to estimate the heat transfer coefficient as a function 
of the pressure ratio and the time, and hence to estimate the heat flow from the 
wall to the gas or vice versa. By an iterative procedure, one could therefore 
obtain a simultaneous solution of the pressure equations (2) and the heat transfer 
equation (44). The appropriate iterative procedure is as follows: 

1. Assume, as a start, a constant gas temperature T g ; calculate the con- 
ductance and outgassing; and upon integration of (2) obtain the pressure 
P', and the density p' from (49) and temperature T g * from (47). 

2. Calculate the coefficient of heat transfer h from (43), and, from (44), 
knowing the wall temperature, evaluate the new gas temperature T . 

K 

3. Calculate as in step one the conductance for t next time interval and 
repeat the procedure. 

The following points concerning the temperature and the isothermal flow 
conditions employed here should be kept in mind: 

1. The temperature is approximately constant during the entire evacuation 
process with the exception of a very brief period at the start of evacua- 
tion. 

2. A complete heat transfer analysis requires the inclusion of all forms of 
heat transfer and of heat flow from all directions. It would have to in- 
clude radiation from space and from other radiating surfaces, conduction 
from adjacent walls, heat sources and other detailed properties of the 
walls. 

3. Techniques for performing complex heat transfer calculations with a 
certain degree of confidence are being devised for spacecraft. Such cal- 
culations could provide temperature versus time for the surfaces which 
could then be used to estimate the gas temperature. 

4. If the temperatures of the walls are known as a function of time, the.r 
effects can be utilized to establish the outgassing rates of surface ma- 
terials and the gas temperature can be estimated and incorporated in 
the gas conductance*. 


METHOD OF SOLUTION - COMPUTER PROGRAM 

The program for the solution of equations (2) was coded in Fortran langu- 
age for the CDC 3100 computer. The numerical integration scheme used was of 
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the fourth order Runge-Kutta type, which requires no starting values other than 
the initial conditions (X Q , Y Q ) and which permits one to change the interval of 
integration, h, at any time without redoing the previous computations. The 
equations used in this procedure are: 


Vi = \ + 1/6 (k 0 + 2k, + 2k 2 + k 3 ) + 0 (h s ) , 


(50) 


where the k are calculated from 


k 0 =hf(X > ,y.) 

k, =hf (X„ + l'2h, Y„ + l/'2k 0 ) 
k, = hf (X n + l'2h, Y n + 1/2 k,) 
k ,=hf (X„ + h, Y„ t k 2 ) 


and 0 (h 5 ) is a truncation error on the order of h s . 


The integration time interval is increased exponentially from its initial 
value, h, to a final value of 2 h in order to speed up the computation and reflect 
the exponential behavior of the pressure time function. To prevent instability of 
the solution, however, the time interval is also modified in the process of inte- 
gration as follows: 

At each integration step, the program evaluates the parameter 


a = 



E 


(51) 


where the k. are the Runge-Kutta variables and E is their average. 
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If a < 1, the time interval used for the integration has provided a value not 
too different from the previous value. The time interval is acceptable and the 
computation continues. 

If cr > 1, the change from one value to the next has been too great and hence 
the interval is unacceptable. As a first step, the interval is reduced by 20%, the 
Runge-Kutta variables are ree validated and the criteria above is rechecked. If 
the criteria that cr be less than one is still unsatisfied, the interval is reduced 
to 1/2 the previous time interval. 

Examples of typical program inputs for the one- and three-volume cases, 
together with the corresponding solutions, arc contained in Appendices A and B 
respectively; the program itself is documented in Appendix C. 


COMPARISON OF EXPERIMENTAL AND COMPUTED RESULTS 

In the process of developing the computer solution to the above system of 
equations with variable coefficients, several combinations of volumes, connecting 
passages, and outgassing materials were investigated. Some of the profiles com- 
puted will be presented in the next section together with an explanation of the 
parameters chosen to characterize these curves. In this section a comparison 
of the calculated profiles with experiment will be made. 

The characteristics of the systems chosen for comparison, such as the 
volume dimensions, the passage dimensions, the type and amount of outgassing 
materials, and the external pressure profile were determined by the experi- 
mental investigation. The tests were also used to establish the pressures at 
.which corona discharge between electrodes in the container would be initiated 
and extinguished and the corresponding times. 

Figure 4 is a schematic of the pressure-time test setup and the orifice 
construction. The one-liter volume was fitted with one of three different ori- 
fices, .437 in diameter, .137 in diameter, and .040 in diameter, connecting the 
compartment to the simulated space environment of the 12-inch diameter by 
18-inch height bell jar. The evacuation was provided by a 140-cfm mechanical 
roughing pump and 2-1500 lit/sec oil-diffusion pumps baffled with liquid nitrogen. 
The bell-jar pressure profile simulates the pressure profile encountered by a 
Thor-Agena B rocket carrying a 1000-pound payload into a 160-mile orbit. A 
Hastings thermocouple pressure gage and a nude ionization gage were used to 
record the pressure in the compartment. The pressure in the bell jar was 
recorded by a Hastings thermocouple gage and an ionization gage. The polymeric 
materials included in the compartment, selected as being typical of the encapsu- 
lants used in flight systems, were: Eccofoam FPH, with a density of 4 lb./cu.ft,; 


33 





12" Dt A. * 18" 


Figure 4-Pressure»Time Test Setup 


silicone rubber, RTV 11; and epoxy, Biggs 823 resin. The test specimens of 
these materials were produced in the laboratory in accordance with accepted 
spacecraft application techniques. The Eccofoam and RTV samples had a sur- 
face area of 140 cm 2 , while the epoxy sample consisted of a conformal coating 
about 0.015 inch thick over an aluminum tube with a surface area of 120 cm 2 . 
The coating was vacuum-processed while curing. The RTV samples were out- 
gassed for 24 hours before the test. Non-outgassed RTV caused fluctuating 
ionization-pressure gage readings in the 10 ~ 4 torr region. This was attributed 
to contamination of the gage elements, which rendered the gage useless for 
pressure measurements. The z4-hdur outgassing ph^se is, however, typical of 
space hardware exposure in environment-chamber tesu. The pressure versus 
time in the bell jar did not follow the same profile in each test but showed a 








variation of about one-half decade because of the influence of the variable gas 
load and orifice sizes. The temperature in these tests was about 21.1°C (70°F). 
Figures 5, 6, and 7 show both experimental and computed results for various 
combinations of orifice size and outgassing material. The pressure profiles for 
each combination were obtained separately and superposed in these figures for 
convenience. 

Comparing these profiles, it is seen that the computed results agree very 
well with the experimental results. The difference (values in the same decade) 
may be attributed to: 

• The difficulty of reproducing exactly the same pressure profile in the 
vacuum jar in the different runs. 

• Gage errors and errors in gage-reading. 

• Differences between actual material properties and those assumed in 
computation. 

For the computation, a mathematical description of the outgassing material is 
used. The u^ta for this representation are obtained from the literature and, 
unless data are available on the specific material being tested, divergence may 
be expected between test and analytical results. 

It is expected that better definition of the materials should produce closer 
results, if needed. Also improvements in computed results can be expected 
when temperature gradients between compartments are included. 


APPLICATIONS AND PARAMETRIC PRESENTATIONS 

Additional pressure profiles for «■ mgle volumes containing various different 
outgassing materials were obtained. These are shown in Figures 8, 9, and 10, 
where they are labeled according to several parametric attributes. Profiles 
for the more complicated three-volume system were also computed, and these 
results are shown in Figure 11, parameterized in the same way. 

The data obtained in the computation of these profiles can be used to esti- 
mate the length of time that should be allowed before certain electrical functions 
which are affected by pressure can be safely performed. They may also be 
used to estimate the size flow passages appropriate, for a given volume or set 
of volumes and given outgassing materials, to insure that suitable pressures are 
reached inside thesevolumee within a specified period of time. It would there- 
fore be very useful to generalize *he results of these computations for application 
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PRESSURE ^TORR) (Log 10 Bos«) 



Figure 5 -Comparison of Experimental and Computed Results for a 1-Liter Compartment 
with Orifice Diameter .437 in., r = 0.1 sec, Ambient temperature (Outgassing material 
types and surface areas are indicated; data on Outgassing rates used in the Computa- 
tion were from ref. 15) 
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Figure 6-Comparison of C/.perimental and Computed Results for a 1-Liter Compartment 
with Orifice Diameter .134 in., r = 1.0 sec, Ambient temperature (Outgassing material 
types and surface areas are indicated; data on Outgassing rates used in the Computa- 
tion were from rjf. 15) 
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Figure 7-Comparison of Experimental and Computed Results fora 1-Liter 
Compartment with Orifice Diameter .040 in., r=10 sec, Ambient Tempera- 
ture (Outgassing material types and surface areas are indicated; data on 
Outgassing rates used in the computation were from ref. 15) 
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to other cases (Ref. 20). This can be done through the use of certain parameters, 
namely those which have been used to label the various curves shown in the figures. 
The choice of these parameters is based on the following considerations. 

The profiles show that within a very short period of time the pressure drops 
to a quasi-steady value. This sharp drop occurs in a region of viscous and 
transitional flow. In this region, the relations expressing the passage con- 
ductance are nonlinear and dependent on pressure and viscosity. Beyond this 
region, the flow is molecular and the conductance depends only on passage 
geometry. The quasi-steady pressure reached and the corresponding evacuation 
time for a given sysrem are thus mainly a function of those parameters which 
characterize the flow in the free molecular regime. In view of this, the profile 
curves have been grouped in accordance with certain parameters in the molecular 
region of the curves which are amenable to generalization, namely the molecular 
time constant, r , the outgassing material gas load at 1 hour vacuum exposure, 

Q, , and the quasi-steady pressure, H . 

The molecular time constant, which relates the size of the volume to the 
pumping speed in the molecular flow region of the passage, represents the time 
needeu for the pressure to drop to l/e = 0.3G8 of the initial pressure and is 


r 



c ) 


(52) 


where V (liters) is the volume of the container to be evacuated and S (Ht/sec) is 
the pumping speed of the passage. The pumping speed for an orifice of area 
A (cm 2) in the molecular flow region (corresponding to a pressure of 10 ~ 2 to 
10 “ 3 torr for the present orifice diameters and for air at 20°C), is a constant: 


S = 1 1 • 6 A (1 /sec ) . 


(53) 


The time constant, which presupposes molecular flow with constant con- 
ductance and no source of gas in the volume being evacuated, is obtained from 
the solution of the linear differential equation expressing the flow balance; i.e., 
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in the molecular region and can be visualized as the slope of the log of the 
pressure profile at time zero. 

If a constant i.e., independent of pressure and time, source of gas (torr- 
1/sec) exists, such as that provided by a leak into the volume, the differential 
equation will be 


2 = sp - 
dt 



(54) 


the solution of which is 


P = P. e" s/v 



(55) 


This equation shows that as t - the first term becomes very small, so that the 
second term controls the pumpdown. For Q L constant, an equilibrium pressure, 



(56) 


is thus obtained. This pressure represents the condition for equilibrium between 
the inflow leakage, Q L , and the outflow through the orifice, SP L . But when the 
in-leakage is caused by material outgassing, the rate of inflow is not constant, 
but is a decreasing function of time. Consequently, after the initial rapid decay, 
the pressure continues to decay slowly, rather than stabilizing at an equilibrium 
value. However, one can define a quasi-steady-state outgassing rate Q L as the 
outflow rate from the material after 1 hour of outgassing in vacuum, and then 
use the corresponding quasi-steady -state equilibrium pressure as one of the 
parameters describing the system. This equilibrium pressure will be: 

= Ql = W (57) 

L 3 S 


% 
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wb^rp (cm 2 ) is the surface area of the outgassing material; 




'tor r 1 i t' 


sec - cm 4 


» 


is the material rate of gas output after 1 hour of vacuum exposure, as given in 
Tables 2 and 3; and S (lit/sec) is the pumping speed for the orifice, as given by 
equation (53). This pressure is not the ultimate pressure obtained in the volume 
but is representative of the slowly decaying pressure. Table 6 lists for various 
systems the geometric data, the outgassing material, the characteristic rate 
Q g , the corresponding load rate Q L , the calculated pumping speed S, and, 
finally, the resulting parameters r and P L . 

CONCLUSIONS 

An analytical tool which provides internal pressure profiles for volumes ex- 
posed to vacuum has been developed. Comparison between calculated and 
experimentally obtained pressure profiles has shown sufficient agreement be- 
tween the two. The computer program developed to provide this information by 
solving a system of differential equations is general and can be applied to a 
system having many volumes, interconnected with various types of passages 
and containing materials outgassing in a known manner, such as is the case in a 
spacecraft. The computer results, obtained in both numer :al and graphical 
form, can be used to provide the following information: 

• The long term pressure conditions in a compartment. 

• The time required to reach the quasi-steady state pressure. 

• The pressure differential across surfaces. 

• The effects of variation of chamber volume and passage size on attain- 
able pressures. 

• The results for an assembly of experiment packages, given the results 
of tests of the individual packages. 

• The equivalent conductance of an experimental package tested in a 
vacuum chamber, including the outgassing effect of materials and sur- 
faces in the package. 
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Table 6 

Volume/Orifice/ Mate rial Combinations and Calculated Parameters 
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•Data from reference <6 



4 The behavior of a system in any other environment, given the results 
obtained in a vacuum chamber. 

Many pressure profiles have been calculated both for a single volume and 
for multi-volume arrangements; some of these profiles have been included 
here. Parameters have been introduced to generalize the results of these 
particular test runs so that they may be used to estimate conditions for other 
similar systems. 


RECOMM ENDATIONS 

Calculation of pressure profiles should be carried out in conjunction with a 
thermal analysis of the system, which provides the transient temperatures of 
the system walls. By the method delineated here, these irrrweratures can bo 
used to estimate the gas temperature and consequently provide 'mproved values 
of the flow passage conductances. The data from the thermal analysis can also 
be used to establish the temperatures of the compartment surface more accu- 
rately. This better knowledge of the outgassing surface temperatures should 
then, in turn, provide a better definition of tho material outgass**ig rates. This 
requires, however, a better knowledge of the variation of mate rial outgassing 
rates with temperature than is available at present. 

As an additional improvement in the calculation of the pressuie profiles, 
the computer program should be made more efficient to provide faster compu- 
tation for very large systems. 
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APPENDIX A 
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TYPICAL COMPUTER INPUT DATA AND CORRESPONDING 
SOLUTION FOR THE SINGLE CHAMBER CASE 
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APPENDIX C 


COMPUTER DOCUMENTATION 


Input Cards 
Program Listing 


Program Flow Charts 
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MAX. ><EC. SIZE 13 **{* 

SUBROUTINE JtUKUTA 

CUNHok I 1HE,T,N,NN*Tm, 1PR|NT*TEmp,qr*ms»H*FTA. ffc-P^PlUJ.XPCWi. 

1 UPtlw)*V<10)«VlNV£N (10 > jNARkJUlfl )*AIIO0A&< 10*5 J « U] UU* > J < UM *0 • > ) » 
1 I AOU 0.5) .nCUNUU.. ) . JC0Nn<lu*S> i iNCfX<10*S> .ui IMI.IU^), 

- * L2lli^iL.5l*t3(lU^llU>) »0UtliSlllL*>n4igiUS2<lXL*:>J*NF'R6;>S' 
i i iM(t>o)» table t5U).NP(ii),AfiEA<iui#is»ch* idelU 

. DiKEASLU* OEU-1A- (2 HA. _ ,_ ... _ . __ 

DATA DELTA / 1 . 0*1 . tJ 1 » 1 . o2,l .03, 1 , "4, 1 , Q5, 1 • U6. 

l i»Q7*l.Le*l .L9*l 7*1.8#:, y »4. 0/ 

DIMENSION Kt<lU)*K2(10)»K3(1U).K4(i0) 

HfcAL R X , K2*J(3,M __ . 

Kt __ 

l* " i * t iht 

H2 « . 

►Ml) = PZER0<T> 

uQ 10A . I « 1>NN . . _ _ . 

1 U *P< i > « p< I ) 

C ALL_ tV AL 

mu 101 1 «. l.N 

i-i mu) » lpuj - . 

f « TiHfc * h2 

. . DQ-10E-i,j -2*AUL 

l 2 XPl U « P( I ) ♦ Kl( 1*1) *H2 

**(1) =■ PZtRUU > _ 

xptl) • P « l ) 

CALL tVAi 

DO 1C.' 1 * 1#N 

iw3 KkU.) * 4)PX4 X _ 

K3 

DO 10“* i * ?*N N 

I- 4 xpti) « pen ^juxx^uam2 

CAUL tVAl 

1)0-105 a- &~1«J L . _ .. 

15 K3t I ) • L>P( I ) 

K 4 

r * i lMt ♦. h 

DO 10b i x 2*NN 

1 6 XP(l) ■ P(l> ♦ JWU_-Xl_*H 

*U> = P2fcRO(T> 

XP(1) « P<1) 

CALL IVAL 

DO 1 u * i = 1* M . _ . 

1 7 K 4 t 1 ) a lip ( [ ) 

rHELK ERROR 
UU 12' i * l.N 

AX MKlti) .♦ K*U)_«> Mm ♦ K4n>) % ,25 

SIU*AbS ( AX -K 1 ( 1 ) > * A0S tAX“K2< I > >*AfcS ( AX *K3 ( I > ) * AH$ <AA»*4U)> 

IF < .«5*bTD7AbS- L* X l_r.- 1-..H X-12 0 ' 1* 0* 13 G 
I'D I'ONTlNOt 

T 1 Mfc = Tint ♦ H .* . 

HO 1 4 u t X ?,NN 

J * l~l - _ . 

AX x 1K1( J)*K2< J)«T2i J)*K3( J) ) * H / 6.Q 
XP( I ) • XP( I ) ♦ AX 


IF IXMI) - XP( 1 ) ) 135# 140*140 

1^0 CON I 1 NUfc- — ... . _ 

DU 145 1 a 2 • NN 

i 5 Pm x xpli^ 

IMintLTA- >d) 141 , 142# 142 

1*1 l DELTA * mill a I . 

1**2 H a M*DfcLTA< IDELl A) 

HfcTufu, — - m 

1 0 IU6LTA a 1 

H x . k*K - _ 

GO TO 1*0 

a. 5 lufcin «- i- 

H a . 5*H 

GO 10 15.0 

END 

it boo fcSCQftfr<S> COP IbO ... 

EOF 
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CALL NUMBER ( 5 , « 1 0 . 1# ♦ 35# FlOATF < JJ ) 

-UEM ? I _ _ 

CALL PLOT(0,#10. #3) 

-Ui). 21 JU X 2-lxUBBMl 
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CALL PLOT<PAF<RAY< 9» I ) » YPAUfc, J Tf N) 
■UEH « St, T-lIfft 


751 CONTINUE 


call plotco. #io..3) 

- UU 732u..I ». 1>1 ABB AT 

Y a PAhRAY ( J, I ) 

u n 

7521 IF CY - l.fc-8) 752o#7522,7522 

-J5Zl .TPAfif >, • t ALflftmXl 


CALL PI 0T(PaRRaY( 9» U*YPAbE.M>> 

_ / 52^ CUNHliyfc 

7504. CALL Pi OT< 10 . . 0 . »*3) 
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i (P( |I«MjH«CU, j) V* lQt| >*E**“T/TAU*QM*ARfcA 


VULUHfc 


COMMON TIMfe.l .N.NN.Tx, JPHM! .TFNP.GRANS.M.Et A. 1 tHP 2 »P( «).XP( 91 . 

1 QPJ 7 ).V(/).* 1 N* ERSE( /).NARfA( 7 J , l 7 # > *»U|I '. 7 ), 

X JCOM^t TT). motxc V>I.CU?.?,*~I.C 20 .'.* ),Ci( 7 , 7,5 >.UM(btt>, 
1 OUTQA$X( 7 » 7 ) » 0 UTQAS 2 ( 7,7 ) . T ABLE 1 50 ) . MP*'ESSU * * fcMHOR. 

1 Oft “A^AU( 7 . 7 V.KCONli< 7 ),Np( ft) • ARt»A < / ) . (SEARCH. I UCI.T A 
IM MERSION 0(1 ») 

HO 1 Q& 1 • ITS 

_ M 1 ) » SOM » 0,0 

API * XP( 1 * 1 ) 

IS. uONTIMUfc 

MA » MARRAM) 


JO IQt M » Ml) 


tMIPRiNT) 101 , 101,110 
l it MPIT6 (61,30 0 ) l,K .QfU,M.TAu(l,lP,Qri I .<*,AHfcAua«> I #K i 
<(k ► OPM A T |5M EVAL.2x,2HI»,i2.3H K«. l2,2X,.»W'i j a ,tX4. 3. * 2H< « , El? . 5. 

1 it. 3 HQF»,fl 2 . 5 . 3 x» 5 H AREA=.fcl?.S) 

10i U( |)s6U) *(OI(1»R)*eKPF(-l/TAUU.ft) >*OM i ,A)l*AHEAO«ASl l#'l) 

1 V_ CUNT ( NUb 

Nu * NCONDt I ) 

U P 102 K « 1 ,MC 


J * jCoNO(i.K) ♦ 1 

1 Q< SUM . SUM ♦ (XP(J)-XPI) ♦ U)KDUCT(I,J,K.XP| .XPI J)> 

upc i ) * ($um*c( m*v inters^ i > 

If < t POINT ) 100 , 1 10,1211 

m m« i tf < 61 . 202 » uoptn 

2 n«. format < 21 h • • • D ERIVATIVE pp< , 12 , 3 h» = ,fel 2 . 5 , 3 * 2 * •) ) 
IK CONTINUE 

UPTURN _ 


END 

FUN CTION CO NDUCT ( ] , J • K , P I , P j ) 

COMMON MMt , T . M.NN.TP, (PRINT , TFNP , BRANS ,H . F T A . 1 ENP 2 ,P< «>,XP< V) , 

Ir OkRiyATlVfe ARRAY PP ( 7 ) HAS SAME NAME AS t«E EOUIVaLENT DIMENSION DP MERE 
1 OPP( 7 ),V« 7 ) *V INVERSE! 7 ),*ARFA< 7 ) , AflFA > 0 *Si 7 # 7 > *Q It 7 , 7 ), 

1 JCO* 0 ±„l>jLL,JjH>B*t J.# 7 ),C'i< 7 , 7 # 5 >,C 2 < 7,>.5 ),CJ< 7 , 7,5 >,TlM( 50 >, 

1 OUTCASH 7 v 7 ), 0 UT 6 A$ 2(/,7 > , TABLE< 90 ) ,NP»ESSUR. ENRON* 

■lj*fX^,2l/T.AUt 7 t 7) # jg;ONjj_/) # MPC ft)#A«6A< 7 ) , I SEARCH, I DELIA 
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C 3 U.J,*) 


P|4 . 4. /J. 1415926536 

V^IXL-U m j-<L^r 2 QD ;.£ 001 v 20®3 > * L 

2u0i 0 * S0RT(P14 * Cl(|» J,X) > 

g»i-T.n. Mte 

2 U(U U * Clt I « J* K ) 

■..y.UL. 10 . 2 fl 06 .., 


2 uqs o * absfick I . J,K ) 


C 2 < 1 . J»K ) ) * .5 


(,0 TO < 1 , 2 , J, 4 , 5 ) ,L 


TYPE 1 * CIRCULAR APERTURE 


1 AP • C] ( !• J,K) 

iLUte-^- « nt i ,mui* urn. 


I0u RATIO * PI/PJ 
U— g-A AIXflL 


lf(AB 5 F(U)-. 97 ) 102 , 102,103 
..XQi-.U -» l./i i 


;o*.u J. 


IMABSFtU) - . 97 ) 102 , 102,104 


-Ut.*.??/ 
J- . 32 ) 


AJALLUK 


102 |F(U- . 32 ) 105 , 106,106 

lO<L_CimOUCT «_23.97»Tft»(U»»*712)«sai| Ta , ?ftft ) * AP/ ( 1 , -y > 
UU TO 1X0 

a.q^- XOO f.lQSf 1PP— 


lov CONDUCT ■ 6 . 27 *TM*AP/tl.-U) 

^ - 14-110 

10 b CONDUCT « 6 . 27 *TM*AP 

. ^ -LJ 1 A/U 2/.112 fc 

111 CONDUCT ■ (CONDUCT ♦ 3 . 62 Mp* Ap )* . 5 

11 *- Mfl TO 1 130 , 1 3 0 , 35 0 ) , L 

lOl CONDUCT « 3 . 62 *TM#AP 


UU 10 ( 150 , 130 , 350 ), L 


15 t UO TO 9999 


TYPE 2 • LQNtt CTRCULAR PIPb 


2 IF (QP • . 015 ) 201 , 200.200 

?Jl. I .t-Uf--, .A» 1 ?Q 2 f 2 lL^ 2 q 3 , 
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K) 
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. SUB ROUTINE QHID . ___ 

l X AND V ARE TM6 PAD! LOCATIONS OP THf LOWER LEFT HAND 

C~” ONE CYUtfe FOR lOQQRIOS, 

l . _ .11 m PR n i s positive, n $ifliari£5 the n. nre.kqf oi v 
c direct ions respectively. 

C ... iF R Oft M IS NEflAl 1 VF < I T _$ IJDttir 1I3LP LQtt. OR1U U 
C AND EQUALS THE NUHBfcft OF CYCLES. 

C . _ Si2£ Of QR1U. EQUALS X$«R |Y YS*N- 
MEAL LOQF 

n .% x 

y2* absf (Floatf(n:> • ys % y 

_ .11 P- Y. . - . 

call plot <x.y»3> 

Jhitti-SJUlSUia ... . 


40 


U i XS 

itul _ 


-ii 


UU 140 I Pl.HM 

_ C*LL.PLQLJ1AjIU*> 

CALL PLOT ( X 1# Y2 . ? ) 

ti i ju r 

BUFF « Yl 



Y2 ■ UUfF 


UO To 160 
At l.:B 


X* a X 

-WLJa* iji,«R 

DO iso J si » 10 

call Plot <xi,vi.2> 
call PLOT ,Xl.Y?.2) 
WIEA-J-Xi. 


Yl 

,11 


Y 2 

BUFF . 


15i X 1 * .4342Q44B*L0QF(FLUATF(J)*1. >* XS ♦ X2 

. X? ,lFLQAtFU))*XS * g 

155 Xl > X? 

i6j; call ..PLOT <ji.v# 3> 
vlsy 

*2 i X 

IMN) 250, *00. 240 
241 1 y2 • YS 
NN*N*1 

U Q-340 _ I «JkNM , 

CALL PLUI ( X i , Y l • 2 ) 

CALL PLOT ( x2 , Yl# 2) 

Yl s Yl* Y2 

k*UFE- .* _Xl_ . _ 

xl « %’i 
34> x2 s BUFF 


UU TO 9 U 0 
25; WB * -N 
Y2 ■ Y 

. - - Y 1_ ■ _ _t . 

UU 355 l »l,N N 



CALL PLOT ( Xl » Yl , 2 ) 
CALL PLOT <X2#Y1.2> 
RUFF ■ xl 


QBIUQ^OU 
CORNER Ct Th,1«IU0u 1U 
» UP THE LENflfl!U0U2^ 
OP1UOU30 
IS l ONS IN IH6QH1U0U40 
QKlUOuSo 
PCSECriYE D I 3« I U 0 U 6 U 
QPIU0v7U 
QKIOOURU 
QHIU0W|> 
QP I UQU90 
QPIU01Q0 
QPJUOUU 
OP1U0129 
QP10013W 
OP I U0l4)i 
QPIU015U 
QPI0016U 
8P1U017U 
QRIOfllftU 
QPIU019U 
QR1002CO 
QP1U021U 
OP1D042U 
QPIU023U 
OPIU024U 
OP1U025U 
0HIU026U 
OP 1 00270 
OHIUO^Ob 
QHI 00*91) 
OPIU040U 
QPJUO^IU 
OHI0032U' 
0HIQQA3U 
3RIU0A4U 
0P1U035U 

OPIU036U 

0H1U047U 

OKI U03bu 

QPl4<U9U 

OP1U040U 
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0P1U042U 

OP 1 1)0430 

OP 10044 

8PIU045U 

0PIJ046U 

QPI0047U 

QP 1 00 480 

OKI 00490 

OPJuoio’iT 

QHIU051U 

QPIU0P2C 

QP 1 1)0)30 

OR 1 00540 
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QRIOO*6U~ 
OPtOOP7u 


x 2 « huff 

J5l Yl « *43429 440 • LOQF < ► LOATf< J> *1 . > • YS »Y2 

Y2 s FLOArF(l)#YS ♦ Y 

3 55 Yl » y2 

9UB RETURN 

END. . ___ _ . . 


.QPimoy __ 

QR 1006 OU 
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QR IDO 630 

QR 1 0*641) 
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mNCT lu N ILOQlO(X> 
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ENTRY mc.tOc 

1 1 H 6 * SUPPORT riil'*TD YIElIj Tfrti OirrERF* 4 XE 
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